arXiv: 1503.05122v2 [hep-ph] 6Jul2015 


Two-loop QCD Corrections to Be Meson Leptonic Decays 


Long-Bin Chen^* and Cong-Feng Qiao^’^^ 

^School of Physics, University of Chinese Academy of Sciences, 
YuQuan Road 19A, Beijing 100049, China 
"^CAS Center for Excellence in Particle Physics, Beijing 100049, China 


Abstract 

The two-loop QCD radiative corrections to the Be meson leptonic decay rate are calculated in 
the framework of NRQCD factorization formalism. Two types of master integrals appearing in 
the calculation are obtained analytically for the first time. We get the short-distance coefficient 
of the leading matrix element to order olf by matching the full perturbative QCD calculation 
results to the corresponding NRQCD ones. The result in this work helps the evaluation of the 
Be leptonic decay constant, as well as the Cabibbo-Kobayashi-Maskawa matrix element \Veb\-, 
to the full next-to-next-to-leading order degree of accuracy. 

PACS number(s): 12.38.Bx, 13.25.Gv, 14.40.Be 


* chenglogbinlO@mails.ucas.ac.cn 
^ qiaocf@ucas.ac.cn, corresponding author 


1 



The advent of non-relativistic Quantum Chromodynamics (NRQCD) factorization 
formalism causes investigations on heavy quarkonium more reliable [l|, which improves 
the understanding of strong interaction. It has been noted that for quarkonium pro¬ 
duction and decays, in many cases the leading order calculation in the framework of 
NRQCD is inadequate. And, mostly the discrepancy between leading order calculation 
and experimental result can be rectihed by including higher order corrections, which has 
stimulated various investigations in this respect. 


Be meson system, which contains two different heavy quark flavors, has some peculiar 
natures different from heavy quarkonium, and recently attracts great interest, especially 
with the progress of the LHCb experiment Q. Though Be meson is very elusive in 
experiment, the feedback from the investigation on it is extremely great, e.g. on some 
aspects of quantum chromodynamics(QCD), weak interaction and even new physics. Of 
the he system, the higher excited states will mostly cascade down to the ground state, 
the pseudoscalar Be meson, through hadronic or electromagnetic transitions, which then 
decays to lighter hadrons or leptons via weak interaction. 

By virtue of nonrelativistic QCD (NRQCD) formalism, the Be meson decay amplitude 
may be expressed as perturbative QCD(pQCD) calculable short-distance coefficients 
multiplied by non-pertubative NRQCD matrix elements. The expression for leading 
order(LO) Be meson leptonic decay width is simple and known for long, and the next- 
to-leading order calculation was completed by Braaten and Fleming two decades ago [^. 
In this work, we compute in pQCD the two-loop radiative corrections to the pseudoscalar 
Be meson leptonic decay rate, i.e., the short-distance coefficient for the leading matrix el¬ 
ement at the next-to-next-to-leading order(NNLO), by matching the perturbative result 
in full QCD with the corresponding perturbative calculation in NRQCD. 


The calculation of massive two-loop Feynman integrals is somehow tough, especially 
with two mass scales. For this reason, there are only a few master integrals with different 
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massive propagators have been accomplished [^-1^. The method of differential equations 
turns out to be an efficient and powerful technique for the calculation of Feynman 
integrals. Recently, it was found by Henn that the solution of differential equation will 
be simplified considerably if the bases of master integrals are chosen properly j^. In 
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our calculation, by employing the technique of differential equation and choosing certain 
bases of master integrals, we successfully obtain the master integrals required in the 
calculation of two-loop QCD corrections to Be leptonic decays. 

The Be meson leptonic decays. Be ^ I I'l with I being e, /i, or r, are heavy-quark- 
annihilation processes through axial-vector current, which are very important to the 
study of Be physics while have not been, but expected to be, measured. Theoretically 
the decay rates can be formulated as: 

, ( 1 ) 

where Vbc denotes the CKM matrix element; M and mi stand for masses of Be meson and 
charged leptons, respectively; and Gi? is the Fermi coupling constant of weak interaction. 
Generally, the Be decay constant is dehned through the transition matrix element 
of charged weak current, as 

(0|67^75c|5e(p)) = ifsy , (2) 


which parameterizes the strong interaction effects and contains both perturbative and 
nonperturbative contributions. 

The short-distance contribution can be isolated and calculated in perturbation the¬ 
ory, by matching the charged weak current in QCD to a series of operators in NRQCD. 
In the rest frame of Be system, up to corrections of order the relative velocity of 
heavy quarks within the meson, the matching relation reads 

(0|67°75c|Sc(p)) = Co {Q\xl^c\Be{p)) + G 2 (0|(DxQt. + ... , (3) 


where Cq and C2 are short-distance coefficients that depend on the heavy quark masses, 
renormalization scale p, and strong interaction coupling The coefficients Co and C2 
will be determined by matching the perturbative calculation of the matrix elements in 
full QCD with what obtained in the framework of NRQCD j^. The coefficient Co was 
obtained at one-loop order in Ref. i.e.. 
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This work is about to calculate the NNLO QCD corrections to the short-distance co¬ 
efficient of the leading order matrix element in expansion analytically, i.e. CO2, by 
which the theoretical prediction for Be leptonic decay rates will come up to the NNLO 
accuracy. 



FIG. 1 : Two loop diagrams that contribute to the concerned processes. For Gi^2,6; the 
symmetric diagrams are implied. The diagram G5 contains contributions from massless 
fermions, gluons, ghosts, and the massive fermions. 


In calculating the two-loop contributions in full QCD, we first revisit the one-loop 
QCD corrections to meson leptonic decays. At one-loop order, there is only one 
diagram in Feynman gauge, and the Mathematica package FeynArts js] is employed 
to generate the amplitude. FeynCalc ^ combined with code written by ourselves are 
used to manipulate the 7-matrix algebra and spin projections, and FIRE together 
with $Apart ll| are employed to reduce all the related integrals into a set of master 
integrals. In one-loop case, the integrals are merely subject to two massive tadpoles. 
After performing the standard renormalization procedure, we then obtain the coefficient 
Co at the order of which agrees with Ref. j3|. 


The topologically independent two-loop order Feynman diagrams are schemetically 
shown in Fig.l, where the solid line represents for bottom quark and dashed line for 
charm quark. Note, within the figure, the solid-dashed lines exchanged diagrams of 
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are implied. Though the boson may couple to charm and bottom quarks 
through both vector and axial-vector currents, in practice only the axial-vector current 
contributes to the annihilation decays of pseudoscalar Be. The calculation of two-loop 
amplitude takes the same procedure as in one-loop case. The hardest part of the two-loop 
calculation in this work resides in the evaluation of the master-integrals, as explained in 
the following. 



fl.3 

FIG. 2 : Two-loop diagrams of various master integrals. 


The master integrals we confront in the calculation are shown in Fig. 2 , where dotted, 
dashed, solid lines correspond respectively to the gluon, charm quark and bottom quark 
propagators, the incoming and outgoing wavy lines of diagrams Is and Jg have a hetitious 
invariant mass of = {mc + mb)‘^. In order to obtain the master integrals by employing 
the technique of differential equation and for the sake of compactness, we attribute the 
two heavy quark masses to a new parameter a; = ^, and normalize the loop integrations 
in terms of 
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with 7 E the Euler constant and /x the renomalization scale. 


The integrals Ji, I 2,13 are just the multiplication of two tadpoles, and the integration 
for tadpole diagram can be found for instance in 12|. The I 4 type of integral can be 
obtained in this procedure: first use AMBRE-package IlSlI to transform the integral 
to Mellin-Barnes representation, then use MB-package [l^ to single out the poles in e 
expansion, and last evaluate the integral by closing a proper contour and summing up 
residues. In the end, we hnd 

^ i ^ ^ ^ ^ ^ ^ ^ 176C(3)) . (6) 

Except for changing mb to rric in the normalization ([7]), the integral I 5 is just the same 
as I 4 , which agrees with Ref. 121. 

Iq can proceed in a similar way as / 4 , while with the coefficient of after the 

reduction. Up to the e^, the Iq reads: 


^ 3 17 .59 7r\ .65 4Qn^ 

.2 

— ( - 16755 + 47507x2 + Mtt^ - 38407x2 ln( 2 ) + 12080C(3)). 
480 ^ 


( 7 ) 


Then integral Jy can be readily obtained the same way as Is. Note that the integral Iq 
will serve as the boundary condition for integrals /io-i 3 - 

The master integrals for /§ to /13 are 
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, ( 11 ) 
( 12 ) 
(13) 


which satisfy differential equations in x. 
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In deriving the differential equations for Ig-n, the FIRE package was employed. We 
hnd the following differential equations exist: 

dis _ W - 10 + {5D - U)x - 6-3D D-2 

dx x{l + x){2 + x) ^~^2 x{2 + x) ^ 2x{l + x){2 + x) 

dh _ 2x{7-3D) D-A + xi2D-5) 2(0-2) 

dx (l + a ;)(2 + a;) ® (l + a ;)(2 + a;) ® a;(l + a ;)(2 + x) ^ 

In Ref. the author suggested that properly chose of master integrals can lead to 

significant simplihcations of the differential equations. To employ this technique, we 
need to find a pair of bases gg and gg, by which the above differential equations can be 
transformed into the following form: 


d^g{e, x) = eA{x)g{e, x) + 0(e, x) . 


(16) 


Here, 0(e, x) is a known function and can be expressed as harmonic polylogarithms 
of argument x. Note, the differential equation (fT 6 H contains poles in x = 0 , — 1 , — 2 . 
Suppose gg = hg{x)Ig + bg{x)Ig, with bg{x) and bg{x) being rational functions of x and 
can be expressed as: 


i=—nQ 

ni 

bgix) = 


rL2 


713 


i=—no 


Vii Vxi Vsi 

X* ^ (1 + x)* ^ (2 + x)* ’ 

2=1 ^ ^ 2=1 ^ ^ 

( 17 ) 

712 713 

Zli 1 Z 2 i 1 Zgi 

'^(1 + 2^)* “^(2 +a:)*’ 

(18) 


by trial and error, we finally sort out two bases that satisfy the canonical form of Eq. llTHll . 
i.e. 


98 = 

99 = 


x^ 

1 + X 


h + h 




2 + 6 x + 5x^ ^ 

2 x ^(1 + x )(2 + x) * 


2 x + 1 

x3(2 + x) ■ 


(19) 

( 20 ) 


Notice the boundary station at x = 1 can be solve the same way as J 4 and Jg, then 
the integrals can be solved iteratively in terms of logarithms and polylogarithms. We 
use HPL-packages 1^ to transform the logarithms and polylogarithms into harmonic 
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polylogarithms with simple argument. In the end, the master integral Ig obtained as 
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With fl2l|) . Jg can be derived out from equation flT^ directly. 


For integrals Jio and In the following differential equations exist: 
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In solving the above differential equations, we determine the bases in a similar way 
as in the case for Jg and Jg, and then the differential equations can also be solved 
iteratively. Note here the logarithmic functions may appear when transforming the 

























differential equations into the form of flT6|) . however the bases can transform the above 
differential equations into a Strictly Triangular Matrix when setting D = A. The bases 
we hnd are 
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Then the analytic form of Jio up to order is obtained by solving the differential 
equations: 
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The Ill hence can be simply obtained from Eq.([22]) and will not be shown here. 


The master integral /12 can be obtained from integral Jio by exchanging rub with rric 
in the norm and replacing x by 1 /x, while the integral /13 achieves in a similar way as 
In. All the analytic results we calculated have been numerically checked with FIESTA 
jl 6 |. The integral 1 12 was given in jd], while we hnd a misprint there that the hrst and 
second terms of (5.20) on page 397 in Ref. jd] should be multiplied by a factor of 1/8, 
otherwise taking the x —)■ 1 limit one can not get the x = 1 result in the same paper. 


After the calculation of master integrals in Fig.l, one can then start the renor- 
malizeation procedure to remove the divergences encountered. The renormalization is 
performed by subtracting the one-loop sub-divergencies and the two-loop overall diver¬ 
gencies. The quark wave functions are renormalized in the on-shell scheme, while the 
strong coupling constant ag is renormalized in the MS scheme. The NNLO renormal¬ 
ization for Be leptonic decays is similar to what shown in Ref. 17|. After performing 
the one-loop mass and coupling constant renormalization, and two-loop wave function 
renormalization, the decay width can be expressed as 




(28) 


which can be expanded perturbatively as: 


T = r“ + aT^' + a^T^' + 0{a^) , (29) 
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ao = a(l + aSZ^l + a‘^5Z‘^l + 0{a^)) . (33) 


Here a and oq denote respectively ^ and is renormalized coupling constant 

and CKso is bare coupling constant. The two-loop renormalized amplitude can be then 
reexpressed as: 
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(34) 
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FIG. 3: Mass-renormalization connter-diagrams. 


For mass renormalization, there are two diagrams shown in Fig.3 to acconnt for 
the one-loop mass renormalization. The one-loop conpling and mass renormalization 
constants as well as two-loop wave fnnctions renormalization constants can be obtained 
from 

After the renormalization procedure, we can classify the two-loop coefficient CO 2 in 
color factors as: 


CO2 — {CpSp-\-C aSa~\~ -\-TfSh)Cf-\-^coi\yl^^ 

4 


(35) 


with 


COI ^ ^ Iiix - 2)Cf 


(36) 


4 ^ 1 -I- x 

being the one-loop coefficient, Np the number of light quarks. The coefficients sp, sa, sp 
and sh read as: 
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(40) 


It is notable that after taking the renomalization procedure in above, the result is 
still divergent, which can be attributed to the anomalous dimension of NRQCD current 
2l| . In MS scheme, the anomalous dirnension of NRQCD current hrst arises at two-loop 

2l| . For our case, the anomalous dimension 


order as noticed and defined in Refs, 
reads: 

d In Zj^jsfjiQCD 


18 


IJ,NRQCD — 


din fi 


1 -I- 6x -1- x2 
2(l + x)2 


2 CfCa-. 2 
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)“s + • (41) 


In the literature, the two-loop QCD corrections to Be leptonic decays, the coefficient 
CO2, was once investigated under the condition of small parameter x = ^ expansion 
up to the second order (^)2 [l8|. For sf,sa and sh we expand our complete analytic 
results in x to the second order and find an agreement with Ref. 18|, while for sl, we 
:ind there is a misprint of redundant term We notice that in the attachment files of 
3 )he renormalization constant Z 2 disagrees with the one given in [3, Hi- have 
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applied our calculation procedure to the study of two-loop corrections to J/'ip and T 
leptonic decays, and found full agreement with Refs, j^, 22|, analytically. 

To evaluate numerically the leptonic decay rates to the next-to-next-to-leading 
order degree of accuracy, we take the following input parameters 


23 


25(: 


rric = 1.5 GeV , rrih = 4.8 GeV , Gp = 1.16637 x 10 ® GeV ^ , 

| 14 fe| = 0.0406 , = 0.499 GeV , = 3 , Aqcd = 0.214 GeV , (42) 

= 0.106 GeV , m^. = 1.777 GeV , t{Bc) = 0.509 ps . 

dere, the non-relativistic decay constant dehned as 


2l| . is obtained by pot ential model evaluation 2^, the pole charm- and bottom-quark 


masses are adopted 25(], and t{Bc) is the B^ life time. For numerical estimates we take 
factorization scale fif = 1 GeV to separate perturbative and nonperturbative domains 
and the renormalization scale is set to be at bottom quark mass. In the end, the decay 
constant in MS factorization scheme can be expressed as: 


= (1 - 1.39( 


C^siX^b) \ no \2\ fNR 


-) -23.7(- 


0 /. 


TT TT 

cNR _ n vnc fNR 


Be 


= (1 - 0.094 - 0.108)/£" = 0.798/£ 


(43) 


With which, we can immediately obtain the branching ratios of B^ leptonic decays at 
NNLO accuracy, i.e.. 


Br{Bc —)■ ~ 1-8 X 10 , Br{Bc —)• /V^z/^) 7.6 x 10 


-5 


(44) 


The branching fraction of Be decay to positron and neutrino is much smaller than the 
numbers in above as expected. 

It is notable that Eq.(jl3]) indicates that the NNLO corrections are greater than the 
NLO corrections, similar as the case of quarkonium leptonic decays. However, the recent 
result on 3-loop corrections for T leptonic decays turns out that the N^LO contribution 
is small and the renormalization scale dependence is greatly reduced 261. 

In summary, we calculated analytically the two-loop QGD corrections to B^ meson 
leptonic decays in the framework of NRQGD. All the master integrals were achieved ana¬ 
lytically by means of Mellin-Barnes integral or Differential Equations. We expanded our 
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analytic results in parameter a; = ^ to the second order and found a partial agreement 
with the results in previous calculation. To conhrm our calculation, we restudied the 
NNLO QCD corrections to heavy quarkonium leptonic decays and can fully reproduce 
those results in the literature. With proper inputs, we numerically computed the Be lep¬ 
tonic decay widths to the full next-to-next-to-leading order degree of accuracy, and found 
the the NNLO corrections are remarkable. This calculation may be helpful as well to 
the precision measurement of Cabibbo-Kobayashi-Maskawa(CKM) matrix element |Wfe| 
when the Be meson decay constant is well determined by lattice QCD calculation. 
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